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Abstract 

It is known that many equations of interest in Mathematical Physics display solutions 
which are only asymptotically invariant under transformations (e.g. scaling and/or 
translations) which are not symmetries of the considered equation. In this note we 
extend the approach to asymptotic symmetries for the analysis of PDEs, to the case 
of difference equations. 

1 Introduction 

When we consider a differential equation the analysis of symmetries leads, via a standard 
procedure [U EH EH HSU E3 > to the determination of invariant solutions. 

This method has been extended to consider conditional symmetries ^3 El HH EHj , 
i.e. transformations which do not leave the equation invariant, but such that some of 
its solutions are left invariant. This notion can be further generalized to consider partial 
symmetries where a subset of solutions, each of them not necessarily invariant, is mapped 
into itself by the transformation at hand jHj (see also [7] for a short review). 

In a recent contribution j^Jj the notion of asymptotic invariant solutions with respect 
to standard symmetries and conditional symmetries for differential equations has been 
introduced. This theory easily extends also to partial symmetries. The approach of [U 
can be seen as a development of |2U| , based in turn on the renormalization group approach 
to differential equations [H EH 03 EH- I n this way we were able to explain the asymptotic 
scaling symmetry observed in numerical experiments analyzing anomalous diffusion and 
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reaction-diffusion equations [321 • I n this note, asymptotic is meant in the sense of "for 
large values of some dependent or independent variables" . 

Another interesting notion in this field is that of approximate symmetries El • Here 
one considers transformations which leave invariant the lower order approximations to the 
solution of a system depending on a small parameter. This notion may also mean an 
approximate symmetry for a series expansion near a fixed point [51 122j. 

All these approaches and results have been introduced in the case of differential equa- 
tions. However, in many physical (or biological, chemical, etc.) models - and in numerical 
simulations of continuous phenomena - one is rather interested in difference equations on 
a lattice, and their symmetry properties ^31 El [Ml El E3 EDI EEj • The reader is referred 
to the bibliography contained in [35] for applications. 

The purpose of the present note is to show that the approach of |21j also extends to 
equations defined on a lattice. 

We will not try to be as general as possible, but rather consider, after a brief intro- 
duction to the generalities of equations on a lattice and their symmetries, a well defined 
equation. We hope that, in this way, the general method will result clearer, and it should 
be easy for the reader to apply it to other equations of interest. In fact, we will focus on 
reaction-diffusion equations 1121 114j , and in particular on the classical FKPP equation 

P3HMHHH]- 

The only possible symmetries for discrete equations if we do not want to change essen- 
tially the underlying lattice are discrete translations, rotations of the lattice variables and 
continuous transformations of the dependent variable. 

We are mainly interested here in equations whose solutions (as observed in numerical 
experiments) are asymptotically well described by a travelling front. We believe that the 
reader interested in more complex asymptotic symmetries will easily understand how to 
deal with them by comparing the present work and [21j . 

In Section [21 we present the generalities on equations defined on a lattice while in 
Section [31 we will consider symmetries for equations defined on a lattice. In Section we 
introduce the FKPP equation both in the continuous and on a lattice and analyze it so as 
to deduce the correct asymptotic behavior by requiring that asymptotically the system has 
the correct required symmetries. Section is devoted to a final discussion of the results 
and some concluding remarks. 



2 Difference equations on a lattice. 

For simplicity we will consider a lattice A in R 2 = (x,t), and a difference equation for 
a real function u(x,t) on it. Higher dimensional lattices and matrix functions defined 
on them could be defined in the same way. We will use the notation of |29[ I38j . Thus 
the points on the lattice, having coordinates (% in ,tm,n)i wm be indexed by a couple of 
numbers (m,n) € Z 2 , and to each site is associated a real variable u m ^ n G R. 

The lattice is intrinsically described by assigning the difference between neighboring 
points and will be of the form 



(1) 
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These will be referred to as the lattice equations. The set of values £, m ,n defines a function 
£ : A — > R via £(xm,n; t m>n , u miU ) = £ TO ,n> an d similarly for 77, r, i?. By the notation u m)Tl we 
mean u(x mj „, i m)n ). The shifts (m,n) —* {m + l,n) and (m,n) — > (m,n + 1) correspond 
to (x — ► x + £(x,i,u),t — > i + r(x,t, u)) and to (x — ► a; + rj(x,t,u),t — » i + $(x,t,u)) 
respectively. The presence of u allows for solution depending lattices. 
We will then consider a difference equation for the variables 

-fm,nM ■ •F\p^m,nt tm,m {^m+j',ri+i}(jj')gZ 2 ] — 0. (2) 

Given some initial boundary conditions, this equation must allow us to get u in all points 
of the lattice (past or future). 

As an example of difference equation we can consider the "discrete heat equation", 
obtained as the natural discretization (but only one of the possible ones) of u t = u xx , 

(St)^ 1 [u m , n+ i - Um^n] = {5xY 2 [u m+ i, n - 2u mi „ + u m _i, n ] . (3) 

where £, m ,n = Sx, rj m;n = 0, T m>n = and $ m)n = St. We can express the discretizations 
in terms of an operator A/%, a difference operator of order j — i. Acting with A/% on an 
arbitrary smooth function fk = F(zk) we have 



3 j 3 




where 6k = Zk+i — Zk- In eq. (J2J) the simplest discrete derivative on t, corresponding to 
j = 1 and i = 0, has been considered. 



3 Symmetries 

A transformation S : (x,t,u(x,t)) — > (x',t',u'(x',t')) will be a symmetry of the lattice 
equation if it takes A C R 2 into itself. It will be a symmetry of the difference equation 
defined on A if it is a symmetry of the lattice and leaves the equation -F m)n [tt] invariant. 
Most of the transformations for the lattice will be discrete symmetries. If we consider 
a uniform regular lattice Ar, the point (m, n) will have coordinates x mn = xo,o + mSx, 
tmn = *o,o + an d the lattice equations are 

Xm+l,n Xm,n — Sx , X m ^ n +i X m n — (5) 
tm+l,n tm,n — , £m,n+l tm,n — St . 

We will denote by p the ratio p = Sx/St (p E R) . On Ajj the discrete symmetries are easy 
to determine. They are given by 

T x : (x, t) (x + Sx, t) [shift in x, equivalent to (m, n) — > (m + 1, n)\; 
Tf : (x, t) — > (x,t + 5t) [shift in t, equivalent to (m, n) — > (m, n + 1)]; 

: (x,t) — > (— x,t) [inversion in x, equivalent to (in, n) — > (— m, n)]; 
i?t : (x,t) — > (x, — t) [inversion in t, equivalent to (m,n) — > (m, — n)]; 
i? : (x,i) — > (—pt,x/p) [rotation by 7r/2 with a scale factor, 

equivalent to (m,n) — > (— n, m)]. 
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Moreover, if we accept SAr C Ar rather than SAr = Ar, we have the further trasforma- 
tion 

S q : (x,t) — > (q\x,q2t) , (91,92) £ Z (discrete scaling) . (6) 

Apart from the discrete symmetries given by the transformation written above we can 
introduce Lie symmetries which are better described in term of an infinitesimal symmetry 
generator X, which, taking into account that the lattice Ar is subject just to discrete 
transformations, will involve only transformations of the dependent variable u m ^ n 

X m ,n — Q{.^m,ni tm,m {^"m+ j ,n+i\ (i j')gZ 2 )^ttm,n ■ CO 

Formula will be the infinitesimal generator of a symmetry for |J2J) if 
1 1 

^ ^ Xm+j-n+jTl^p = (8) 
j=-l i=0 

Eq. (JSJ) is equivalent to the request that the flow generated by eq. Q, when u m ^ n = 
Um,nW, is compatible with 

^m,n,\ — Q{%rn,m tm.ni {^m+J,ra+i}(ij)eZ 2 ) (9) 

Let us notice that the symmetries generated by the infinitesimal generator 0, if ^ 
(0,0), are not Lie point symmetries but generalized symmetries. 

Solutions invariant with respect to the symmetries of infinitesimal generator (|2J) are 
obtained by solving the difference equation (J2J) together with the invariance condition 
Q{x 

m.m tm,ni {^ i m+j,n+i}(ij')gz 2 ) — 0. A particularly interesting class of function symme- 
tries is when the function Q takes the form 

which in the continuous limit goes over to point symmetries. 

When the transformation is not a symmetry, but there is an invariant solution, we say 
that we have a conditional symmetry for our problem. In this case to get the invariant 
solution we can add to the equation the condition Q = 0. 

If we consider an equation -F mi n[u] which does not explicitly depend on x and i, it is 
invariant under T x and Tt transformation. In this case we can consider travelling wave 
solutions, depending on Q = x — vt. Defining v = k(5x/5t) with k € Z, the travelling wave 
solutions can be obtained by imposing [23| 



Q = (f- k - f t )u m , n . (11) 

For k = 1 this can be read, similarly as in the continuous case where the condition is 
u x + (l/v)ut = 0, as Q = A~u m ^ n + ^/Afu m ^ n where A~ and are the 'down' and 'up' 
discrete derivatives defined as 



(12) 
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As the regular lattice An is invariant under 

M = f£f t : (m,n,u miTl ) -> (m + k,n + l, it m +fc,n+l) , 



we can pass to the variables (( = x — vt,r = t,w((,r) = u(x — vt,t)). Obviously Q and 
M are related, as we have Q = Tt{M~ x - 1) u mn . So we can introduce two new indices 
(fj, = m — kn, v = n). In these new independent variables the difference equation F mtn [u] 
will transform into a new difference equation for the variable w^^, ^^^[it^j. 

The reduction of eq. (j2J by the symmetry (fTTj) implies that w^ )U+ i = w^y, so, as w^ )U 
is independent on u, we can write w^ v = and reduce ^^[ui] to $^[^4]. If this equation 
admits a solution, it represents a travelling wave solution. 

Let us now consider a continuous transformation on A generated by the infinitesimal 
symmetry operator 

X = *A^, (13) 

where ^ is a constant. The infinitesimal transformation (|13j) corresponds to the transfor- 
mation A^ = A^ + \$>A^ + 0(A 2 ). If the quantity 

W» := (M^A]) - 

is not zero but goes to zero for large fj,, then X is an asymptotic symmetry of 3>„L4]. 

Finally, if there exist A-invariant asymptotic travelling wave solutions of the original 
equation F min [u], we will say that X is Bin Gsyuiptotic coTiditioncil syTYiTuctvy of -F^^^u]. 



4 The discrete FKPP equation and its symmetries 

In the rest of this note we will focus on the standard FKPP |15| I25[ I33| reaction-diffusion 
equation 

u t = u xx + u(l - u). (14) 

See e.g. El for a discussion of its properties. Here we are specially interested in its 

asymptotic solutions for large t. For sufficiently localized non-negative initial data, 

x->_+oc -ko(x+ki) (h. u \ > n for r > 
u(x,t = 0) = { ' ^0' fel ^ >u ' Iora; > u (15) 

^ =1, for x < 0. 

the asymptotic solutions correspond to a stable front travelling with constant speed and 
smoothly connecting the "fresh" (unstable) stationary state u = ahead of it and the 
"exhaust" (stable) stationary state u = 1 behind it. 

The difference equation representing its simplest discretization on a regular lattice A# 
is given by 

F m ,n{u\ : ■ j t = (16) 
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Let us stress that Sx and St are small but finite parameters, representing the lattice spacing 
in the x- and in the t-direction respectively. 

As noted above, in analyzing the symmetries of eq. , we should restrict our attention 
to transformations S : (x,t,u(x,t)) — > (x',t',u'(x',t')) which are also symmetries of the 
lattice equations ©• 

Let us first consider transformations which do not act on the u m)n . It is easy to see 
that shifts in x and t are symmetries of the equation (|16|) , This is also the case for 
the inversion in x, but not for the inversion in t. Also, since x and t are intrinsically 
different here, the rotation R is not a symmetry of the equation. We can recover the time 
inversion by considering the map B t : (x,t,u(x,t)) — > (x,—t,u(x,—t)), corresponding to 
(m,n,u mtn ) — > (m, —n,u m - n ). This represents a time inversion and backward dynamics. 

It is easy to see (by the same argument as for the continuous equation) that there is 
no way to have a scaling symmetry of the independent variables for eq. ()16j) . 

4.1 The discrete FKPP in a moving frame 

Following the results presented in Section 13 a set of values u m ^ n will represent a travelling 
solution with speed v = k5x/5t = kp (with k G Z) if, for any j G Z, 

^m+kj,n+j — H"m,n • (1^) 

The equivalent of passing to a moving frame of reference will be passing to the independent 
variables 

C = x - vt ; t = t. (18) 
The inverse of the change (|18|) is given by t = r, x = ( + vt. It follows from 

Cm,n = x m n -vt m n = m(5x) - ^-n5t = (m - kn)5x , 

ot 

T~rn,n — tm,n — Tl dt . (19) 

that the lattice equations, in terms of the new coordinates, are given by 

Cm+l,n Cm,n — ^X , Cm,n+1 Cm,n — V dt — k bx , 
1~m+l,n T"m,n — , Tm,n+1 T"m,n — St . (20) 

The lattice point with coordinates (£ = ^Sx,t = vSt), has coordinates (x = mSx,t = 
nSt) in the original system, with fj, = m — kn, v = n. The inverse map is m = /j, + kv, 
n = v. 

We will denote by w^^ the dependent variable associated to the point of coordinates 
C = fJ.Sx, t = uSt. We have 

We can now use (|2T|) to express (fTB|) in the new variable w. We get 

w m—k(n+l),n+l ~ w m—kn,n 



— W m — kn,n{\ W m —kn,n) ~\~ (^2) 
+ 



St 

tVm+l—kn,n ^^m—kn.n ^m—l—kn,n 

Sx 2 
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Passing now to the indices (/i, v), we can rewrite eq. (|22|1 as 

: — ^ = — g£ ~ + w^ u (l-w^ u ). (23) 

Eq. (|23|) can be simplified by recalling that 5t = -Sx. So the equation Q^ v \w\ reads 

[wn+i,v ~ 2iVix,v + «V-i,"] ~ Z Sx l w n-k,u+i ~ v>n,v] + Sx 2 w fl)U (1 - w^ u ) . (24) 

Let us recall that 5x is a small but finite parameter, representing the lattice spacing in 
the x-direction. 

If we reduce eq. 1)24(1 with respect to a proper combination of the translations, we 
get that the solution of 1(24(1 will depend only on the £ (not on the r) coordinate, i.e. 
w fi,v+j = w fM,v f° r an J € Z. If it^j, does not depend on the r variable, we can as well 
restrict to a fixed value of the v index, and write simply A, := Wn U . Inserting this ansatz 
into 1)24(1 we get the reduced difference equation 



[A t ,+i-2A IM + A^ 1 ] + -fe[^-^_ fc ] + 8x 2 A lx {l-A IJ ) = 0. (25) 
4.2 Travelling fronts 

A solution to 1(25(1 is, by construction, a travelling wave solution to the difference FKPP 
equation 1)16(1 . In the continuum case, we know this will actually be a travelling front. In 
the moving frame coordinates (£, r), it will quickly tend to zero for positive £ and to one for 
negative Q. Moreover, it is stable. If we start with a sufficiently regular initial datum, the 
(non translation-invariant) solution will tend to this moving front as time goes by. Finally, 
the front solution is characterized by a sharp transition between the two stationary states 
u = and u = 1. The region with small u is described by an exponential, u ~ exp (— aC) 
(with 1/aa characteristic width). 

The equivalent of the exponential behaviour in our discrete setting is 

A^+t = (l-a5x)A IM . (26) 

In fact, as fi = Q/5x, linife^o A^ = j4oexp (—«£). Solving eq. (|2l)|) we get 

Ap = (1 — a5xYA . (27) 

Putting formula 1)27)1 in 1)25)1 . with fc = 1, we have 

5x 2 [a 2 - (u + &c)a + !-(!- afe)^ +1 A ] = 0. (28) 



Eq. 1)28)1 is not solvable for fi — > oo if a < 0. The case a > gives us a monotonically 
decreasing solution for growing fi. For large /i, we can leave out the last term of eq. (jSBj) 
and solve it for a: 

a 2 - (y + fe)a + 1 = 0. (29) 

This equation admits two real solutions for v + Sx > 2; studying the stability of the 
equation as in the continuous case, one finds that the velocity is v = 2 — 6x. So the front, 
observed in the numerical integration of the discrete FKPP on Ar, would be just a little 
slower (of the order of 5x) than the front, which is solution of the FKPP equation. 
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Note that we have just found an asymptotic symmetry. In fact the choice of a scaling 
invariant solution (f2l)|) is equivalent to set 

Q = aA fl + A C A M (30) 

The corresponding transformations, a scaling on A and a translation on £, are not a 
symmetry of the non- linear equation (|25|) as = X^^A] — ^^[XA] = — a<5x 3 (l — 
a5x)Aj l+1 ^ 0. As for a < A^ is a decreasing function of fi, we get that lim^oo = 0. 

5 Discussion and Conclusions 

In this note we have applied the method of asymptotic symmetries to a discrete version of 
the FPPK equation. In such a way we have shown that also in the discrete case the FKPP 
equation possesses an asymptotic scaling symmetry which gives the correct asymptotic 
behaviour of the system. 

Long time ago Ibragimov and Fushchich ^ ^] introduced the notion of approximate 
symmetry (for differential equations); this is based on some small parameter e appearing 
in the equation, and is a symmetry "up to terms of order e". One could wonder what 
is the relation between approximate symmetries in this sense and our approach to the 
front solutions of the FKPP equation; note however that no explicit small parameter e 
is present in the FKPP equation. A natural way to introduce such an e would be to 
look for solutions whose amplitude is of order e, i.e. u„ jm = ev n ^ m (with |« n ,m| of order 
one); proceeding in this way, however, the obtained scaling invariant solution will be 
small everywhere, not only asymptotically. Actually, in our approach the natural small 
parameter is inversely proportional to the group parameter generated by the asymptotic 
symmetry. Thus the small parameter comes into play only when we consider an additional 
object (the asymptotic symmetry). It appears, therefore, that our notion of asymptotic 
symmetry cannot be reconducted to the notion of approximate symmetries of Ibragimov 
and Fushchich ^ ^]; further work would be needed for a better understanding of the 
relation between these two notions. 

Work is in progress to obtain the asymptotic solutions for the anomalous reaction 
diffusion equation. This equation is nonautonomous in the independent variables, so that 
no point symmetry is available. In such a case the use of a renormalization of the lattice, 
provided by eq. ©, might be important. 
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